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Abstract 

Based on a criterion of mathematical simplicity and consistency 
with empirical market data, a stochastic volatility model has been ob- 
tained with the volatility process driven by fractional noise. Depend- 
ing on whether the stochasticity generators of log-price and volatility 
are independent or are the same, two versions of the model are ob- 
tained with different leverage behavior. Here, the no-arbitrage and 
completeness properties of the models are studied. 
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1 Introduction 

In liquid markets the autocorrelation of price changes decays to negligible 
values in a few minutes, consistent with the absence of long term statistical 
arbitrage. Because innovations of a martingale are uncorrelated, there is a 
strong suggestion that it is a process of this type that should be used to 
model the stochastic part of the returns process. As a consequence, classical 
Mathematical Finance has, for a long time, been based on the assumption 
that the price process of market securities may be approximated by geometric 
Brownian motion 

dS t = fiStdt + aS t dB (t) (1) 
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Geometric Brownian motion (GBM) models the absence of linear correla- 
tions, but otherwise has some serious shortcomings. It does not reproduce 
the empirical leptokurtosis nor does it explain why nonlinear functions of 
the returns exhibit significant positive autocorrelation. For example, there 
is volatility clustering, with large returns expected to be followed by large 
returns and small returns by small returns (of either sign). This, together 
with the fact that autocorrelations of volatility measures decline very slowly 
DP) 0) has the clear implication that long memory effects should some- 
how be represented in the process and this is not included in the geometric 
Brownian motion hypothesis. The existence of an essential memory compo- 
nent is also clear from the failure of reconstruction of a Gibbs measure and 
the need to use chains with complete connections in the phenomenological 
reconstruction of the market process DQ. 

As pointed out by Engle [5], when the future is uncertain investors are less 
likely to invest. Therefore uncertainty (volatility) would have to be changing 
over time. The conclusion is that a dynamical model for volatility is needed 
and a in Eq.flT]), rather than being a constant, becomes itself a process. This 
idea led to many deterministic and stochastic models for the volatility ([0], 
[7] and references therein). 

The stochastic volatility models that were proposed described some par- 
tial features of the market data. For example leptokurtosis is easy to fit but 
the long memory effects are much harder. On the other hand, and in contrast 
with GBM, some of the phenomenological fittings of historical volatility lack 
the kind of nice mathematical properties needed to develop the tools of math- 
ematical finance. In an attempt to obtain a model that is both consistent 
with the data and mathematically sound, a new approach was developed in 
[S]. Starting only with some criteria of mathematical simplicity, the basic 
idea was to let the data itself tell us what the processes should be. 

The basic hypothesis for the model construction were: 

(HI) The log-price process log St belongs to a probability product space 
(fli xf^2, P1XP2) of which the (f2i, Pi) is the Wiener space and the second one, 
(f^2) P2), is a probability space to be reconstructed from the data. Denote 
by oj\ G Vii and u 2 G fi 2 the elements (sample paths) in Qi and f2 2 and by 
T\ )t and Tif, the cr-algebras in Qi and f2 2 generated by the processes up to 
t. Then, a particular realization of the log-price process is denoted 

logS 4 (wi,w 2 ) 

This first hypothesis is really not limitative. Even if none of the non-trivial 
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stochastic features of the log-price were to be captured by Brownian motion, 
that would simply mean that St was a trivial function in f^. 

(H2) The second hypothesis is stronger, although natural. It is assumed 
that for each fixed u>2, log St (-,^2) is a square integrable random variable in 

These principles and a careful analysis of the market data led, in an 
essentially unique wa}0, to the following model: 

dS t = fi t S t dt + a t S t dB(t) (2) 

\oga t = /3 + - {B H (t) — Bjj (t — 5)} (3) 


the data suggesting values of H in the range 0.8 — 0.9. In this coupled 
stochastic system, in addition to a mean value, volatility is driven by frac- 
tional noise. Notice that this empirically based model is different from the 
usual stochastic volatility models which assume the volatility to follow an 
arithmetic or geometric Brownian process. Also in Comte and Renault [9] 
and Hu [10J, it is fractional Brownian motion that drives the volatility, not 
its derivative (fractional noise). 5 is the observation scale of the process. In 
the 6 — > limit the driving process would be a distribution-valued process. 

The equation (j3J) leads to 

at = ee l{B H (t)-B H (t-8)}-^) 2 8^ ( 4) 

with E [a t ] = 6>0. 

The model has been shown [8] to describe well the statistics of price re- 
turns for a large 5-range and a new option pricing formula, with "smile" 
deviations from Black-Scholes, was also obtained. An agent-based interpre- 
tation pj] also led to the conclusion that the statistics generated by the 
model was consistent with the limit order book price setting mechanism. 

To our surprise, this data-reconstructed model was met with some hostil- 
ity by the Mathematical Finance community. Perhaps in part because much 
of the nice results in this field are based on simple GBM and also because 
fractional Brownian has been associated to arbitrage. In fact, in the past, 
several authors tried to describe the long memory effect by replacing in the 
price process Brownian motion by fractional Brownian motion with H > 1/2. 

1 Essentially unique in the sense that the empiricaly reconstructed volatility process is 
the simplest one, consistent with the scaling properties of the data. 
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However it was soon realized [12], [32], [33], [15] that this replacement im- 
plied the existence of arbitrage. These results might be avoided either by 
restricting the class of trading strategies [16J, introducing transaction costs 
[T7j or replacing pathwise integration by a different type of integration [18] 
[19]. However this is not free of problems because the Skorohod integral 
approach requires the use of a Wick product either on the portfolio or on 
the self-financing condition, leading to unreasonable situations from the eco- 
nomic point of view (for example positive portfolio with negative Wick value, 
etc.) [20]. 

The fractional volatility model in Eqs.(j2Jl3]) is not affected by these con- 
siderations, because it is the volatility process that is driven by fractional 
noise, not the price process. In fact a no-arbitrage result may be proven. 
This is no surprise because our requirement (H2) that, for each sample path 
oj 2 £ ^2, log St (-, u) 2 ) is a square integrable random variable in fii already 
implies that J a t dB t is a martingale. The square integrability is also essential 
to guarantee the possibility of reconstruction of the a process from the data, 
because it implies [2TJ 

f (^j = Ht{-,u 2 )dt + a t {-,oj 2 )dB t (5) 

Our aim in this paper is to give a solid mathematical construction of 
the fractional volatility model, discussing existence questions, arbitrage and 
market completeness. 



2 No-arbitrage and incompleteness 

Let (fix, F\, Pi) be the complete filtered Wiener probability space, carrying a 
standard Brownian motion B = (B t ) Q<t<oo and a filtration Fi = (^ r i,i) <t<oo- 
Let also (f2 2 , J-" 2 , P 2 ) be another probability space associated to a fracFional 
Brownian motion Bh with Hurst parameter H £ (0, 1) and a filtration F 2 = 
o<t<oo g enera t e d by B H . 
Let us now embed these two probability spaces in a product space (fl, J 7 , P) , 
where Q is the Cartesian product Qi x Q 2 and P is the product measure 
Pi <S> P 2 - We also introduce 7Ti and 7r 2 , the projections of Q onto Qi and 
Q 2 , as well as the a— algebra M generated by all null sets from the product 
cr— algebra T\ ® J- 2 , that is, 

N = a ({F C ^ x n 2 \3G £ J"i ® T 2 such that F C G and (Pi <g> P 2 ) (G) = 0}) . 
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Moreover, we let T = (T\ eg) T 2 ) V Af, the a— algebra generated by the union 
of the a— algebras T\ <g> To and A/". Then F = (Tt)^ is the filtration for 

° \ / 0<t<oo 

T t = {T 1; t®T 2>t )vM. _ 

Furthermore, we extend B and Bh to F— adapted processes on (Q, T, ~P) 
by B 0J2) = (B o 7Ti) (ui,u 2 ) and B H (ui^) —_{B H o 7r 2 ) (u x , u 2 ) for 
(wi,^) G Q. Then, it is easy to prove that B and Bh are Brownian and 
fractional Brownian motions with respect to P and are independent. For 
notational simplicity, hereafter B and Bh will stand for B and Bh- 

We now consider a market with a risk-free asset with dynamics given by 



dA t = rA t dt A = 1 



(6) 



with r > constant and a risky asset with price process S t given by Eqs.([2])- 
OH]), with /x t a F-adapted process with continuous paths, k the volatility 
intensity parameter and 6 the observation time scale of the process. 
Eq.© leads to Eq.gD 



with 6* = Ep [<r t ] , hence 07 is a measurable and an F— adapted process satis- 
fying for all < t < 00 



E 7 
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f{B H {s)-B H {s-5)} 



ds 



using Fubini's theorem and the moment generating function of the Gaussian 
random variable B H (s) — B H (s — 8). 

Moreover f Q \fi s \ds being finite P— almost surely for < t < 00, an 
application of Ito's formula yields 



S t = Soexp 



fi s - 7} a s^j ds + J a s dB s 



Lemma 2.1. Consider the measurable, ¥ -adapted process defined by 

r - {it 



It 



< t < 00 



Ot 



(7) 
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with /i G L°° ([0,T] x flj and denote by r\ = (77 
nential of ( J* ^ s dB, 



c 



the stochastic expo- 



0<*<oo' 
f 



that is, 



Then, 



7] t = exp 



1 



'jsdBs / %ds \ , < t < 00 



1 r T 

o / ls 
z Jo 



lids 



< 00, < T < 00 



(8) 



(9) 



Proof: We make of use the fact that the sample paths of the fractional 
Brownian motion Bh are Holder continuous of any order a > strictly less 
than H, that is, there is C a > such that, P— almost surely, \Bh {t) — 
Bh (s) I < C a \t — s\ a , for every t, s G [0, 00). Then, for all < T < 00 

cT 



E 7 



< E 7 



2 Jo 



exp 



exp 



l s ds 



n k 2 5 2H ~ 2 rT 
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(r + \\iJt, s 



2# 2 



2 r T 

2 ^ 2H - 2 ' e -2|(B( S )-B( S -«5))^ t 



loo; c k 2 5 







< exp < T(r+ 2 ^' L) e^ 2 "' 2 -^- 1 > < 00 



Additionally, we assume that investors are allowed to trade only up to 
some fixed finite planning horizon T > 0. 

Proposition 2.2. The market defined by |D] ; (TJ|) and (Ej) is free of arbitrage 
Proof: Because, by the lemma above the process 7 in ([7]) satisfies the 
Novikov condition (Q, the nonnegative continuous supermartingale rj in ([8]) 
is a true P— martingale. Hence we can define for each < T < 00 a new 
probability measure Qt on Tt by 

—— = Vt , P- a.s. (10) 

Then, by the Cameron-Martin-Girsanov theorem, for each fixed T G [0, 00), 
the process 



B* = B t - 



T —^ds 0<t<T 

°s 



(11) 
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is a Brownian motion on the probability space (fl, J-'t, Qt) ■ 
Consider now the discounted price process 

Zt=~r 0<t<T 
At 

Under the new probability measure Qt, equivalent to P on Tt, its dynamics 
is given by 

° S Z S dB* s (12) 

o 

and is a martingale in the probability space (TI^t^Qt) with respect to 
the filtration [J't) Q<t<T - ^ ^ ne fundamental theorem of asset pricing, the 
existence of an equivalent martingale measure for Z t implies that there are 
no arbitrages, that is, Eq t [ZtlJ 7 ,,] = Z s for < s < t < T. ■ 
Another important concept is market completeness. We note that, in 
this financial model, trading takes place only in the stock and in the money 
market and, as a consequence, volatility risk cannot be hedged. Hence, since 
there are more sources of risk than tradable assets, in this model, the market 
is incomplete, as proved in the next proposition. 

Proposition 2.3. The market defined by (TJP ; ([3|) and (Ti|) is incomplete 

Proof: Here we use an integral representation for the fractional Brownian 
motion [22], [23J 

B H (t)= [ K H (t,s)dW s (13) 
Jo 

Wt being a Brownian motion with respect to P, independent from B t and K 
is the square integrable kernel 

K H (t, s) = C H s^- H I (u - s) H ~^u H -^ du, s <t 

J s 

(H > 1/2). Then the process 



V ' t = exp \W t - -tj (14) 

is a square-integrable P— martingale. Then, defining a standard bi- dimensional 
Brownian motion, 

w; = (B t ,w t ) 
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the process rf t = r\ t r\ t 



r? t = ex P [J v s • dw* - i y ||r s || 2 ds 

where, by lemma 2.1, V = (7,1) satisfies the Novikov condition, is also a 
P— martingale. Then, by the Cameron-Martin-Girsanov theorem, the pro- 
cess 

w; = (w* (1 \w* {2) 



defined by 



W t * {1) = B t - I ls ds 







w t * {2) = w t 



is a bi- dimensional Brownian motion on the probability space (fl, Tt, Qt) > 
where Qt is the probability measure 

= Vt = VtVt (15) 

Moreover, the discounted price process Z remains a martingale with respect 
to the new measure Q^. being an equivalent martingale measure distinct 
from Qt, the market is incomplete. 

■ 

As stated above, incompleteness of the market is a reflection of the fact 
that there are two different sources of risk and only one of the risks is tradable. 
A choice of measure is how one fixes the volatility risk premium. 



3 Leverage and completeness 

The following nonlinear correlation of the returns 

L (r) = (\r(t + r)\ 2 r (*)) - (\r (t + r)| 2 > (r (t)> (16) 

is called leverage and the leverage effect is the fact that, for r > 0, L (r) 
starts from a negative value whose modulus decays to zero whereas for r < 
it has almost negligible values. In the form of Eqs. §Z§ the volatility process 
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a t affects the log-price, but is not affected by it. Therefore, in its simplest 
form the fractional volatility model contains no leverage effect. 

Leverage may, however, be implemented in the model in a simple way [21] 
if one identifies the Brownian processes B t and Wt in (T5]) and ([TBI . Identifying 
the random generator of the log-price process with the stochastic integrator 
of the volatility, at least a part of the leverage effect is taken into account. 

The identification of the two Brownian processes means that now, instead 
of two, there is only one source of risk. Hence it is probable that in this case 
completeness of the market might be achieved. 

Identifying the random generators, a new fractional volatility model is 
defined [23] 

dS t = HtS t dt + a t S t dW t 

rt 3 

log<r t = + k' {t-uf-*dW u (17) 

J — oo 

Proposition 3.1. The market defined by [Tty and (E)) is free of arbitrage 
and complete. 

Proof: The proof of the first part of the proposition is analogous to that 
of proposition 2.2. In fact, a similar argument to Lemma 2.2 yields that 

is a Px-martingale with respect to (J~i,t) < t< T and the probability measure 
Qt, defined by = ^ T 1S an equivalent martingale measure. 

Now that we have shown that the set of equivalent local martingale mea- 
sures for the market is non-empty, let Q* be an element in this set. Then, 
recalling that (^ r i,t) <i<T ^ s * ne augmentation of the natural filtration of the 
Brownian motion W t , by the Girsanov converse [25] [26] there is a {Fi,t) a<t<T 
progressively measurable M-valued process such that the Radon-Nikodym 
density of Q* with respect to P\ equals 

Moreover the process W£ given by 

w; = w t - [ <p s ds 

Jo 
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is a standard Q*— Brownian motion and the discounted price process Z sat- 
isfies the following stochastic differential equation 



dZ t = (/it - r + a t (p t )Z t dt + a t Z t dW^ 

Because Z t is a Q*— martingale, then it must be hold /i(t, u)—r+a(t, u)<f)(t, uj) 
almost everywhere w.r.t. meas x P in [0,T] x f2, meas being the Lebesgue 
measure on the line. It implies 

n T-n(t,u) 

(j){t,U) 



a(t, uj) 



a. e. (t,co) G [0, T] x f^. Hence Q* = Q, that is, Q is the unique equivalent 
martingale measure. This market model is complete. 



4 Remarks and conclusions 

1) Partially reconstructed from empirical data, the fractional volatility model 
describes well the statistics of returns. The fact that, once the parameters 
are adjusted by the data for a particular observation time scale S, it describes 
well different time lags is related to the fact that the volatility is driven not 
by fractional Brownian motion but its increments. 

Specific trader strategies and psychology should play a role on market 
crisis and bubbles. However, the fact that in the fractional volatility model 
the same set of parameters would describe very different markets [8] seems to 
imply that the market statistical behavior (in normal days) is more influenced 
by the nature of the financial institutions (the double auction process) than 
by the traders strategies [IT] . Therefore some kind of universality of the 
statistical behavior of the bulk data throughout different markets would not 
be surprising. 

The identification of the Brownian process of the log-price with the one 
that generates the fractional noise driving the volatility, introduces an asym- 
metric coupling between at and St that is also exhibited by the market data. 

2) In this paper, mathematical consistency of, both versions, of the frac- 
tional volatility model has been established. This and its better consistency 
with the experimental data, makes it a candidate to replace geometrical 
Brownian as the standard market model. 
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